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Dominant Multipoles in WMAP5 Mosaic Data 
Correlation Maps 
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The method of correlation mapping on the full sphere is used to study the properties of the ILC 
map, as well as the dust and synchrotron background components. An anomalous correlation of 
the components with the ILC map in the main plane and in the poles of the ecliptic and equatorial 
coordinate systems was discovered. Apart from the bias, a dominant quadrupole contribution in 
the power spectrum of the mosaic correlation maps was found in the pixel correlation histogram. 
Various causes of the anomalous signal are discussed. 
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1. INTRODUCTION 

Non-Gaussian features for the low harmonics, the 
quadrupole and octupole, were found almost im- 
mediately after the first WMAP mission data 
release [UEllS]. These features are expressed in the 
planarity and alignment of the two harmonics, which 
was noted by Tegmark et al. 4^ . The axis along which 
the extrema of these multipoles are lined up in the 
WMAP mission map was later named the Axis of 
Evil [5]. Note that this axis is not the only selected 
axis in the mission data, the naturally selected axis 
lies in the Galactic plane. In the further releases of 
the WMAP mission data [71 [3 [9] the existence of 
the Axis of Evil was also confirmed. 

Different estimations of the significance of exis- 
tence of this axis, and several hypotheses on its origin 
were made. Various studies, e.g. [HlllO] investigated 
the contribution of background components and their 
influence on the alignment of multipoles {£ = 2 and 
£ — 3), and indicated a small probability of the back- 
ground effect on the orientation of the low multipoles. 
In [11] , where the multipole vectors were used for the 
estimates of this effect, it was also noted that the posi- 
tions of the quadrupole and octupole axes correspond 
to the geometry and direction of motion of the Solar 
System and are perpendicular to the ecliptic plane 
and the plane, given by the direction to the dipole. 
Randomness of such an effect is estimated by the au- 
thors as unlikely at the signiflcance level exceeding 
98% and exclude the effect of residual contribution 
of background components. Continuing the research 



done, Copi et al. [12] conclude that the characteristics 
of low multipoles are abnormally different from ran- 
dom, which may be due to the statistical anisotropy of 
the universe at large scales, or to the problems of the 
ILC (Internal Linear Combination) signal deconvolu- 
tion method. Park et al. [13] note that the planarity 
of the quadrupole and octupole is not statistically 
significant. They also stress that the residual photon 
radiation in the ILC map does not affect significantly 
the level of the effect. 

Cosmological models were developed to explain 
the prominence of the axis in the orientation of mul- 
tipoles. The alignment of the quadrupole and oc- 
tupole could be explained within the framework of 
these models. Various models include the anisotropic 
expansion of the Universe, rotation and magnetic 
field [HITiKTl]. 

Other studies investigated the contribution of 
ecliptic dust in the microwave background. Diego et 
al. [17] have constructed a combination of initial data 
of the WMAP channels: V + W - 2Q, which was 
smoothed over a 7-degree diagram. The CMB is not 
included in this combination, however, it can be used 
as a tool to investigate the behavior of the large-scale 
residual signal. The authors have shown the presence 
of radiation in this residual signal. Its distribution is 
correlated with the position of the ecliptic plane. The 
paper investigates the contribution of zodiacal dust 
and shows that the radiation level of this contribu- 
tion is not sufficient to explain the level of residual 
signal in the map combinations. In addition, a contri- 
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bution from unresolved sources was estimated, which 
could also become an explanation of the observed ef- 
fect. The authors do not consider this contribution 
sufRcient as well. 

Dikarev et al. [18] used the WMAP data and the 
infrared survey data to compute the parameters of 
dust particles, which could provide the correspond- 
ing contribution in the microwave background. The 
authors believe that the silicate beads a few millime- 
ters in size in the zodiacal cloud and trans-Neptunian 
region with optical depth on the order of 10~^ could 
explain additional radiation at the level of 10 /iK in 
the microwave background in the ecliptic. 

We investigated the problem of correlation prop- 
erties of the WMAP maps, given the fact [12] that in- 
creased correlations were previously reported in var- 
ious ILC signal distribution regions. In addition, it 
was proved [201 1211 122] that a great contribution to 
the non-Gaussianity of data on various multipoles is 
given by the background components. Coming back 
to [19], let us note that the distribution of the corre- 
lation coefficients in mosaic areas of the map had a 
great shift for the correlations with the dust compo- 
nent of the background. This may be due to a more 
complex model of the dust distribution than the one 
currently used for the separation of components, and 
then the effect may occur in the correlation at other 
galactic latitudes. 

In this paper, we continue to study the detected 
correlation [19] at different angular scales. To study 
the effect, we calculate the angular power spectrum 
of the mosaic correlation map and select its main har- 
monic components. We apply this analysis to the ILC 
map correlations both with the dust component, and 
with the synchrotron radiation. 



2. MOSAIC 
METHOD 



CORRELATION 



To analyze the map's properties on different angular 
scales, we expand the signal distributed on the sphere 
into the spherical harmonics (multipoles): 



(1) 



= 1 m=-e 



where AS* {9, (j)) are the signal variations on the sphere 
in polar coordinates, i is the multipole number, m is 
the multipole mode number. The spherical harmonics 
are determined as 



where x = cos 0, and P"'(a;) are the associated Legen- 
dre polynomials. For a continuous function A5'(x, 4>) 



the expansion coefficients aim are expressed as 

pl n2TT 

aim = / dx / A5(x, (/))Yf*j(x, 0)ci0, (3) 



where denotes a complex conjugation Yim- 

The correlation properties of two maps of the sig- 
nal distribution on the sphere can be described, on 
a given angular scale, by a correlation coefficient for 
the corresponding multipole i as: 

K{i) = - 7=-' ^ , 

m——£ m——£ 

where tgrn and sirn are the variations of two signals in 
a harmonic representation. The coefficient K{i) can 
be used to assess the correlation between the harmon- 
ics on the sphere, i.e., to compare the properties of 
maps on a given angular scale. However, in the case 
of a search for correlated areas, which do not repeat 
in other regions of the sphere, this approach smears 
such single areas in the process of averaging over the 
sphere within a certain harmonic. In this case it be- 
comes practically impossible to identify the correlated 
areas. 

Verkhodanov et al. [19] proposed an approach, 
which was implemented in the second release of the 
GLESP package 53j . The proposed procedure makes 
it possible to find correlations between two maps in 
the areas of a certain angular size. In this method, 
each pixel with number p subtending the solid an- 
gle Sp is assigned the cross-correlation coefficient 
K (Eplimax) between the data of the two maps on the 
corresponding area. Thus a correlation map is con- 
structed for two signals T and S, where the value 
of each pixel p [p — l,2,...,iVo, and A^'o is the total 
number of pixels on the sphere) with the subtending 
angle Sp , and computed for the sphere maps with the 
initial resolution determined by imax is equal to 



^ ^_ (T(e.,0,)-r(2p))(S(e„0,)-5(Hp)) 



(4) 



Here T(0i, (pj) is the value of the signal T in the pixel 
with the coordinates {9i,(j)j) for the initial resolution 
of the pixclization of the sphere; S{6i, (f>j) is the value 
of the other signal in the same area; T(Sp) and >S'(Sp) 
are the mean values averaged over the area Sp and 
obtained from the higher resolution map data deter- 
mined by imax, fTp, and as^, are the corresponding 
standard deviations in the area considered. 
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Figure 1: Maps of correlation coefficients on the 
sphere, constructed for mosaic correlations of the ILC 
and dust maps. Top to bottom: maps for the pixeliza- 
tion with the pixel side of 16ff , SOff and 54(f ■ 

3. ILC AND DUST BACKGROUND 
COMPONENT MAP CORRELA- 
TION 

In [in] we found a significant bias in the distribution 
of the correlation coefficients for the WMAP ILC and 
dust maps. Let us view this distribution in detail. To 
do this, we construct mosaic correlation maps with 
pixels having the sides of 160, 300 and 540 minutes 
of arc with the corresponding maximal multipoles 
imax — 33, 17 and 9 in accordance with the definition 
(4). The resulting maps are shown in Fig. 1. 

The distribution of pixels (correlation coefficients) 
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Figure 2: Histograms of correlation coefficient distri- 
bution for the ILC and dust maps. Top to bottom: 
histograms for the pixelization with the pixel side of 
160, 300 and 54ff . The dotted line, short and long 
dashes mark the zbfj, ±2f7 and :h3a levels in the pixel 
distribution for correlating random signals, computed 
for the ACDM cosmological model and dust. 
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Figure 3: Power spectra C{£) of the correlation coef- 
ficient maps of the ILC signals and dust. Top to bot- 
tom: power spectra for the pixelization with the pixel 
side of 160, 300 and 54^ff . The dotted line, short and 
long dashes mark the iter, ±2ct and ±3a levels in the 
scatter of signal values for the correlation of a random 
map and dust. 



for the maps in Fig. 1 is demonstrated in Fig. 2. Here 
the ±cr, ±2cr and ±3(7 levels in the pixel distribu- 
tion for the correlation of a random signal and dust 
are marked with dotted, short and long dashed lines, 
respectively. These levels were computed by averag- 



ing the results of calculations for the 200 random 
Gaussian maps modelled in terms of the ACDM cos- 
mology power spectrum. Apart from the shift in the 
distribution to the negative values, there is a signif- 
icant distortion of the distribution shape, which dis- 
tinguishes it from a normal distribution. The median 
values of the distributions for the correlation pixels 
of 160', 300' and 540' are equal to -0.2187, -0.2326 
and —0.2736, respectively. It should be noted though 
that the resulting distribution bias is similar to the 
bias, obtained at the statistical evaluations of the 
quadrupole signal restoration quality using the ILC 
method [24j . 

The shift of the distribution in the negative direc- 
tion defines the monopole in the correlation map. It 
is most likely related to the procedure of component 
separation and is also observed in the correlations of 
one-dimensional cross sections [25l[26]. However, the 
distortion may be caused by the presence of enhanced 
correlations of higher multipoles in the maps. To in- 
vestigate this phenomenon let us construct the angu- 
lar power spectrum: 



Ci£) 



1 



\a.eo\ 



(6) 



The power spectra of the maps with the corre- 
sponding imax — 33, 17 and 9 are demonstrated in 
Fig. 3. The dotted, short and long dashed lines here 
mark the cr, 2a and 3(T levels in the scatter of spec- 
trum values for the correlation of random maps and 
dust. The monopole signal has previously been re- 
moved. It is clear from Fig. 3 that for all the three 
power spectra the quadrupole is much stronger than 
the other harmonics. Using the expansion by the for- 
mula (1), we extract this harmonic from the correla- 
tion map data. The results are demonstrated in Fig. 4. 

A peculiarity of this quadrupole is its sensitivity to 
both coordinate grids we selected: ecliptic and equa- 
torial, as demonstrated in Fig. 5. The figure shows 
that the cold spots near the galactic poles lie in the 
ecliptic plane, and hot spots lie exactly in the equa- 
torial poles. 

Note as well the fourth harmonic in the power 
spectrum for these correlations with the scale of 160'. 
Its amplitude is also higher than the 3cr level. It is 
evident from the map of this multipole (Fig. 6) that 
this correlation is determined by the structures in the 
Galactic plane, which, in general, is anticipated on the 
scales of 15°-20°. 

Returning to the quadrupole, we show that the 
removal of the Galactic plane from the calculations 
does not affect the position of this harmonic. For this 
purpose, we use the KpO mask, which is available from 
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Figure 5: Ecliptic (left) and equatorial (right) coordinate grids in the quadrupole map, constructed based on 
the data with the pixelization of 16ff . The coordinate grid poles fall on the hot spots (light gray), and the cold 
spots (dark gray) lie in the equator. 





Figure 6: Map of the fourth multipole with a superim- 
posed galactic coordinate grid. The map is calculated 
for the correlation window with the resolution of Idff . 





Figure 4: The quadrupole extracted from the ILC and 
dust correlation maps. Top to bottom: maps for the 
pixelization of correlations within the areas with the 
sides of 160, 300 and 54ff . 



the WMAP mission websit^ and is applied for shield- 
ing the Galactic plane and point source radiation. 

Figure? shows the map of dust and ILC signal 
correlations in the areas with a side of 160' with a 
superimposed KpO mask. It also shows a histogram of 
counts outside the mask. Although a full power spec- 
trum has no physical meaning in this case, we still use 
it to assess the statistical properties of the signal, ap- 
plying a set of simulated maps with a superimposed 
mask. A statistical analysis can be then done based 
on the comparative characteristics. Figure 7 shows the 
angular power spectrum calculated on the full sphere, 
where the pixels located in the KpO mask area have 
the zero values inscribed. The dotted and dashed lines 
show the a, 2a and 3cr levels in the scatter of the sig- 
nal values for the correlation of a random map and 
dust with a superimposed mask. Note that the ratio 
of the quadrupole amplitude in the spectrum in the 
studied correlation map to the 3tT level is close to that 
anticipated from random maps, without the use of the 
mask it is equal to 2.2, while in the case of mask this 
ratio amounts to 5.5. Therefore shielding the Galac- 
tic plane increases the significance of the quadrupole 



^ |http : / /lambda . gsf c . nasa ■ gov | 
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Figure 7: ILC and dust signal correlation map for the 16C/ correlation window and its statistical diagrams. 
The KpO mask is superimposed on the map. The diagrams were calculated accounting for the mask. Top left: 
a correlation map for the pixel side of 16Cf . Top right: a histogram of pixel distribution in this map outside 
the KpO mask. The dotted, short and long dashed lines mark the a, 2a and 3a levels in the scatter of the 
signal values for the correlation of a random map and dust with a superimposed mask. Bottom left: the power 
spectrum. Dotted and dashed lines demonstrate the scatter, estimated similarly to the previous figure. Bottom 
right: the quadrupole map. 



harmonic contribution in the correlation map. Fig- 
ure? as weh demonstrates the quadrupole map, the 
axis of maxima of which is turned with respect to 
that shown in Fig. 4 due to the presence of the mask 
in the spectral analysis, and the axis of minima is 
unchanged. 

4. CORRELATIONS BETWEEN ILC 
AND SYNCHROTRON BACK- 
GROUND COMPONENT MAPS 

In [19] we discussed the ILC and synchrotron ra- 
diation signal correlation maps. These maps are as 
well visually revealing the Galactic plane (Fig. 8), but 
lack a clear cut deviation of the histogram from the 
Gaussian distribution (Fig. 9), as it is the case of the 
ILG and dust emission signal distribution correlations 
(Fig. 2). But note that for all the three count distri- 
bution histograms, constructed for the 160, 300 and 
540 'windows, there still exists a deviation from the 
Gaussian case above the la level. 

Let us build the power spectra of the mosaic cor- 
relation maps for the ILC signal and synchrotron ra- 
diation in the same manner as we previously did for 
the dust maps. They are demonstrated in Fig. 10 for 
various correlation windows. Each figure also shows 
the limits of the signal scatter at the icr, ±2cr and 
±3(7 levels, calculated from the data of 200 correlation 



maps for the random signal in the ACDM cosmology. 
We pre-subtracted the monopole from the maps with 
the 160 and 300'correlation windows. 

The spectral analysis shows that similarly to the 
case of dust, there is a very notable harmonic i = 2. 
In order to study its orientation let us build its map 
(Fig. 11). 

Figure 11 draws our attention to the location of 
a harmonic, the orientation of which does not co- 
incide with the Galactic plane. Although hot spots 
are considerably shifted relative to the maxima of the 
correlation maps with the dust (Fig. 4) , the positions 
of cold spots (minima) indicate a notable sensitivity 
both to the ecliptic plane (Fig. 11, left, a quadrupole 
for the scale of 160'), and to the equatorial plane 
(Fig. 11, right, a quadrupole for the scale of 540'). 
Both minima in the two maps lie at the correspond- 
ing zero latitude. 

5. DISCUSSION 

A stable presence of the quadrupole axes linking the 
hot and cold spots, located in the same areas of the 
sky for different ILC and dust component correla- 
tion map windows in conjunction with the power 
spectrum, exceeding \a manyfold, indicates a high 
level of non-Gaussianity of low multipoles in the 
WMAP signal. A distinctive feature of this paper 
from those of other authors as well demonstrating 
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Figure 8: Maps of correlation coefficients on the 
sphere, constructed for mosaic correlations of the ILC 
and synchrotron radiation maps. Top to bottom: maps 
for the pixelization with the pixel side of 16ff , 300 and 
540'. 



the non-Gaussianity of the WMAP low muhipoles 
(see Introduction) is not only the method employed, 
but two more components: (1) we show that when 
cleaning the WMAP data applying the ILC tech- 
nique, the dust component gives a strong anticorre- 
lation to the emitted CMB, manifested both in the 
distribution of correlation coefficients, and in the an- 
gular power spectrum. The synchrotron component 
is considerably manifested within this approach only 
in the power spectrum; (2) the distribution of cor- 
relation coefficients allows speculating about a sort 
of a special signal that has the "knowledge" of not 



only the ecliptic coordinate system, but even more so 
about the equatorial grid. 

Note that two negative spots detected in the cor- 
relation map's quadrupole are composed of negative 
values, indicating the inverse behavior of the selected 
ILC map and background components in these areas. 
The negative areas of the quadrupole strictly corre- 
spond to the position of positive residual signal zones 
detected in [10] . This fact confirms the prominence 
of these WMAP map areas. An important argument 
from our point of view is the position of these spots, 
which does not seem to be entirely random. The coor- 
dinates of the minima of the quadrupole's cold spots 
are close to the (0,0) and (180°, 0) coordinates both in 
the ecliptic, and in the equatorial coordinate systems, 
and the positions of hot spots (the poles either in the 
ecliptic, or equatorial systems) are sensitive to the 
size of the selected correlation scale. We have demon- 
strated this effect plotting a grid on the quadrupole 
map (Fig. 5 and 11). 

Note that we can not yet confidently tell what has 
caused the sensitivity to the coordinate system. We 
can enumerate several potential effects that may play 
a role in enhancing the signal in the ecliptic: 

• the presence of the dust component 
in the ecliptic plane; despite the weakness 
of contribution of the zodiacal dust |10], we 
can not as yet rule out the effect of silicone 
condensates [TS] : 

• the presence of the microwave background, de- 
termined by the solar wind; the results of research in 
this direction, presented in the form of signal distribu- 
tion maps on the full sphere, are not yet available; in 
addition, such a background has to be variable owing 
to the solar activity; 

• the residual noise associated with the inhomo- 
geneous sensitivity of a WMAP pixel, depending on 
the distance from the ecliptic plane; the produced mo- 
saic correlation and modeling for the WMAP noise 
map has shown very low correlation characteristics 
with the ILC map (Fig. 12); 

• The influence of huge dust belts of giant plan- 
ets on the effects of component separation of the 
microwave radiation, the remnants of which are no- 
ticeable in the ecliptic plane. The giant dust ring, 
newly discovered by NASA's Spitzer infrared space 
telescop^ [37], with the size larger than one degree 
should give an unaccounted contribution to the radia- 
tion in the ecliptic plane as a result of multiple passes 
through the power beam of the WMAP antenna. The 
level of such a contribution in the microwave range 



^ http : //science .n asa.gov/headllnes/y2009/ I 
/07oct_giantring . htm 
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Figure 12: Left: pixel histogram for the mosaic correlation map of ILC and noise. Right: power spectrum of 
this map. The correlations were computed inside the areas with the pixel side of 16C/ . The dotted line, short 
and long dashes mark the ±a, ±2cr and ±3a levels in the pixel distribution for correlating random and noise 
signals. 



was not evaluated. 

A discovery of the equatorial coordinate system 
in the correlations, apparent in the positions of the 
quadrupole hot and cold spots, remains a surprise. 
The possibility of its physical manifestation in the 
Lagrangian point L2 seems doubtful. There is no 
evidence of the existence of the Earth's giant dust 
ring, similar to those discovered around Saturn, as it 
would be detectable both from Earth and from satel- 
lites. Neither do we find any signatures of a giant 
tail of the Earth's magnetosphere stretching up to 
L2, which could host areas of microwave radiation. In 
this case, the effect of systematics, possibly associated 
with the process of observations is probable, but not 
clear. If we assume that non-Gaussian deviations of 
the low multipole characteristics can be related to the 
problem of the component separation, including the 
poorly-studied radiation of the ecliptic plane, then 
the question of how these correlations turn out to 
be sensitive to the equatorial system requires further 
study. 

The data anticipated from the Planck mission is 
of much interest in this regard. These data, owing to 
higher resolution and sensitivity, will allow to com- 
pute and investigate mosaic correlation maps in the 
ecliptic and equatorial planes. 
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Figure 9: Histograms of correlation coefficient distri- 
bution for the ILC and synchrotron radiation maps. 
Top to bottom: histograms for the pixelization with 
the pixel side of 16Cf , 30ff and 54-(f . The dotted line, 
short and long dashes mark the icr, ±2cr and ±3a 
levels in the pixel distribution for correlating random 
signals and synchrotron radiation. The random signal 
is computed for the A CDM cosmological model. 
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Figure 10: Power spectra C{£) of the correlation coef- 
ficient maps of the ILC signals and synchrotron radia- 
tion. Top to bottom: power spectra for the pixelization 
with the pixel side of 16Cf , 3 Off and 54ff ■ The dotted 
line, short and long dashes mark the ±a, ±2a and 
±3a levels in the scatter of signal values for the cor- 
relation of a random map and synchrotron radiation. 



Figure 11: The quadrupole, extracted from the ILC 
and synchrotron radiation correlation maps. Top to 
bottom: maps for the pixelization of correlations 
within the areas luith the sides of 16ff , 30 ff and BJ^ff . 
The first map has the ecliptic coordinate grid super- 
imposed, the second and the third — equatorial. The 
coordinale grid poles fall on the hot spots (light gray), 
and the cold spots (dark gray) lie in the equator. 



